Phase transitions of bipartite entanglement 
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We study a random matrix model for the statistical properties of the purity of a bipartite quantum 
system at a finite (fictitious) temperature. This enables us to write the generating function for the 
cumulants, for both balanced and unbalanced bipartitions. It also unveils an unexpected feature 
of the system, namely the existence of two phase transitions, characterized by different spectra of 
the density matrices. One of the critical phases is described by the statistical mechanics of random 
surfaces, the other is a second-order phase transition. 
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The bipartite entanglement of small quantum systems 
(such as a pair of qubits) can be given a quantitative 
characterization in terms of several physically equivalent 
measures, such as entropy and concurrence PQ. The prob- 
lem becomes more complicated for larger systems and/or 
higher dimensional qudits [2] . The interest of character- 
izing entanglement for these systems is twofold: on one 
hand, it has fascinating links with complexity [3] and a 
related definition of multipartite entanglement [J ; on the 
other hand, it has applications in quantum information 
and related fields of investigation [5] . 

In this Letter we intend to characterize the statistics 
of the entanglement of a large quantum system. We shall 
tackle this problem by studying a random matrix model 
that describes the statistical properties of the purity of 
a bipartite quantum system. In the context of quan- 
tum information this model was introduced in [6j [8] in 
order to describe the statistics of the eigenvalues of the 
reduced density matrix of a subsystem and extract the 
first moments of some quantities of interest , like the en- 
tanglement entropy or the purity. We will obtain the 
exact generating function of the purity in the limit of 
large space dimension (large N in the matrix model) and 
will connect the entropy with the volume of the manifolds 
with constant purity (iso-purity manifolds) . We will also 
show that the matrix model undergoes two phase transi- 
tions, one at a negative and one at a positive (fictitious) 
temperature. The phase transition at negative tempera- 
ture will be paralleled to another one, that is well known 
in the study of random matrix models and conformal 
field theory literature [5]. We notice that techniques re- 
lated to those presented in this Letter have been recently 
employed [10] to analyze the statistics of the lowest eigen- 
value of the reduced density matrix. 

Consider a bipartite system in the Hilbert space TL = 



Ha ® U B , with dim Ha = N < dimH B = M. Assume 
that the system is in a pure state |?/>) G H. The reduced 
density matrix of subsystem A reads 

PA = Tr B \iP)(ip\ (1) 

and is a hermitian, positive, unit-trace N x N matrix. 
Its purity 

kab = Ti~aPa e [1/N, 1] (2) 

is a good measure of the entanglement between the two 
subsystems: its minimum is attained when all the eigen- 
values are — 1/N (completely mixed state, maximal en- 
tanglement between the two bipartitions), while its max- 
imum detects a factorized state (no entanglement). We 
consider a typical pure state \ip) [3|8], sampled accord- 
ing to the unique, unitarily invariant Haar measure. The 
significance of this measure can be understood in the fol- 
lowing way: consider a state vector |^o) an d let consider 
a unitary transformation — U\ipo)- I n the least set 
of assumptions on U, the measure can be chosen ran- 
domly in a unique way. The final state l^) will hence 
be distributed according to the Haar measure mentioned 
above (indipendently of IV'o))- Notice the analogy with 
the maximum entropy argument in classical statistical 
mechanics. By tracing over subsystem J5, this measure 
translates into the measure over the space of Hermitian, 
positive matrices of unit trace [5J [5] 

dp(p A ) = V P A(detpA) M - N S(l-Trp A ), 

= d N xj[(x i -x j ) 2 Y[xf8(i-j2^),m 

i<j £ k 

where Xt are the positive eigenvalues of pa (Schmidt co- 
efficients), we dropped the volume of the SU(N) group 
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(which is irrelevant for our purposes) and pN = M — N 
is the difference between the dimensions of the Hilbert 
spaces Ha and TLb- 

We will consider the statistical properties of the 
rescaled quantity 



R = Rab = N 3 tt A b- 



(4) 



The moments of this function can be obtained by lengthy, 
direct calculations. We will propose a different approach 
that makes use of a partition function: 



Z A b = J dp{p A ) exp {~PR A b) 



(5) 



where [3 is a fictitious temperature. This approach is 
easily generalizable to any other measure of entangle- 
ment. The fictitious temperature in the partition func- 
tion (which is the generating function of the purity) is a 
"tool" to fix the value of the purity (and thus of entangle- 
ment). In particular for (3 = one obtains typical states, 
while for larger values of /3 one gets more entangled states 
(for /3 — > oo maximally entangled states) . 

Henceforth we willassume N ^> 1. We will analyze 
in detail the case p = and then give the results for 
M — N = pN > 0. Our problem has been translated 
into the study of random (reduced) density matrices pA 
with 

Zab= f d N XT[(X i - A,) 2 <5(l-f>)e-^ 3 M. 

(6) 

As a first step, we introduce a Lagrange multiplier for 
the delta function 



Zai 



.V- I ^ f d N X 

2n Jx,>o 



xe 



(7) 



By assuming N large we can look for the stationary point 
of the exponent with respect to both the A^'s and £. The 
contour of integration for £ lies on the real axis but we will 
soon see that the saddle point for £ lies on the imaginary 
£ axis. It is then understood that the contour needs to 
be deformed to pass by this point parallel to the line of 
steepest descent. The saddle point equations are 

-2/3iV 3 A, + 2V-^— -z7V 2 e-0, (8) 
' V " Xj 

£Ai = l- (9) 



In the limit of large N, by adopting the natural scaling 
A, - jj\(xi), < x t = ^ < 1, (10) 



we can write Eq. (pi as 



(3X- 



v , P(A') 



where 



p(X) = / dxd(X-X(x)) 
Jo 



(11) 



(12) 



is the density of eigenvalues. A similar equation, re- 
stricted at (3 = 0, was studied by Page [5]. 

We start at high temperatures ^ < 1 and assume a 
solution of the form [T3] (see Fig. [IJ 



a — X 



A 



(13) 



for < A < a and otherwise. This form satisfies the 
integral equation as can be promptly verified. The La- 
grange multiplier £ is related to the parameters a, b by 
£ = i(3{a — 6), and it is purely imaginary, as anticipated. 

We can find a, b by imposing normalization and the 
constraint, which derive from (12 1 and ([9|, 



dXp(X) = 1, 



f dXp(X)X=l. 
Jo 



By imposing the form (13 1 we find 

P 



P ,2 



1G 



; a(a + 26) = 1, 
For /?_ < (3 < (3 + with 

/3_ = -2/27, /3+ 



a 2 (a + 6) = 1. 



(14) 



(15) 



(16) 



there is a unique solution of these equations that yields 
real, positive p(X) : 



(17) 



13a 2 

where A = (^-/3//3_ + ^1 - Z?//?,) 1 / 3 . Notice that 

a(/3)~ 4-8/3, 6(/3) - /T 1 -4/3, for /3 -> 0, 
a(/3)~2 + 6(/3), &(/?) ~ (/3+ - /3)/4, for /3 T /3+ 
a(/3) ~ 18 + &(£), tor pip-, 

HP) 



-12- Vl2(l-/3//3_), 
The average purity is given by 



for/3j/3_. (18) 



N 3 



(19) 



By using (l8| one shows that i?(/3 = 0) = 27V 2 , R(P+) = 
5N 2 /4 and R(P-) = 9N 2 /A (see later for the significance 
of this values) . 



One can also compute the free energy 

F = R — / dx dy\o E \X{x) - X{y)\ (20) 

P Jo Jo 



and using the saddle point equations (11) it is possible 
to show that 



d\p(X) / dA'p(A')log|A' - A| 



= / ,/A/XAi [logX + p~ + i^X) . (.211 



where we also used (14 1, and obtain 



F 

iV2 



1 



(6 — a)a — 



2 + alog(a/4) 3a 4 /3 



a/3 



256 



(22) 



in terms of the function a((3) introduced above. 

Notice that (3F is the generating function for the con- 
nected correlations of R. The radius of convergence in 
the expansion around (3 = defines the behavior of the 
late terms in the correlations. 

One can find the values of all the cumulants of R, ttab 
(or connected correlations, the derivatives of logics) 
in the unbiased distribution at f3 = 0, when p(X) = 
(l/2-7r)y(4 — A)/A. One starts by observing that a se- 
ries expansion of (Tl7|) yields 



«(/?) = £ 4^3 



1-32 



(31-1)! 



/>() 



(21 + 1)1(1 -1)! V/3 







(23) 



By making use of this expression one finds 



(-1)™ d n 



N 3n dp 



(PF) 



2 n+1 (3n-3)! 

N 3n-2 (2n)! 

(24) 

The first three cumulants are of course the large- N limits 
of known results [6 (for small N exact expressions for the 
first 5 cumulants can be also found in [7]). 

We are now ready to unveil the presence of two phase 
transitions. The most evident one is at the end of the 
radius of convergence of the small (3 expansion, which 
occurrs at /?_. We can extend our equations smoothly 
down to (3- but not below. At /3_ we have p(X) = 
2/(27tt)(6 - A) 3 / 2 /%/A and ir AB = 9/AN (see Figures [lj 
and [2]) . The derivative at the right edge of eigenvalue 
density vanishes and some eigenvalues can evaporate to 
+oo.|15j The limits (3 — > /3_ and N — > 00 can be com- 
bined (double-scaling limit) to interpret the free energy 
as the partition function of random 2-D surfaces (a the- 
ory of pure gravity) . Using (24 1 we see that around (3- 
the free energy F oc ((3 — (3Jyt 2 + less singular [9]. In 
fact, if one relaxes the unit trace condition, our partition 
function Z has been studied in the context of random 
matrix theories |llj . The objects generated in this way 
correspond to chequered polygonations of surfaces. Our 




FIG. 1: Density of eigenvalues at different temperatures. The 
phase transitions occur at (3+ = 2 and at j3- = —2/27. 



calculations show that the constraint Tr p A — 1 is irrel- 
evant for the critical exponents. 

The other phase transition occurs as (3 is increased (the 
temperature decreased). The value of b decreases con- 
tinuously and eventually vanishes at /3+ (where ttab — 
5/ AN), becoming b < for (3 > /?+. The solution ([l3l 



is not valid anymore, since p(X) becomes negative for 
A < —6/2. We have to look for another solution, and, by 
noting that at + , p(X) = (/3 + /tt)^A(2 - A) (see Fig. [[J, 
we do so in the usual semicircle form 



p(X) = -VA - bVa - X. 



The normalization and the constraint yield 



(25) 



{3 



(a-b) 2 = l, 
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{a-bf{a + b) = 1. (26) 



This can be easily solved to find 




and hence 



R = N 2 1 



1 



Moreover, from ( 20 1-( 21 ) , one gets 



F 

N 2 



= 1 



3 

AP 



1 

2ft 



log(2/3). 



(27) 



(28) 



(29) 



We can now notice how the phase transition at j3+ is due 
to the restoration of a Z2 symmetry P ('parity') present 
in Eq. (Ill, namely the reflection of the distribution p(X) 
around the center of its support (A = a/2 for (3 < (3 + and 
1 for (3 > f3 + ). For (3 < (3 + there are two solutions linked 
by this symmetry, and we picked the one with the lowest 
F; at /3+ this two solutions coincide with the semicircle 
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( 25 1 , which is invariant under P and becomes the valid 



and stable solution for higher (3. 

One can also determine the expression for the entropy 
S — (3(R — F), which counts the number of states with 
a given value of the purity. The expression for f3 < /3+ 
is quite involved and we will not write it here, while for 
(3 > (3+ it is easy to see that: 



In the critical region, (3 — > /3 + , we find 



TV 2 



log 2 



6((3 - (3+) 



16 



(30) 



,(31) 



where 9 is the step function. We see that S is continuous 
at the phase transitions, together with its first derivative 
although the second derivative is discontinuous. So this 
is a second order phase transition. 

Notice that the entropy is unbounded from below when 
(3 — > +oo. The interpretation of this result is quite 
straightforward: the minimum value of ttab is reached 
on a sub-manifold (isomorphic to SU(N)/Zm [TJ]) of di- 
mension jV 2 — 1, as opposed to the typical case vectors 
which form a manifold of dimension 2N 2 — N — 1 in the 
Hilbert space Tl. Since this manifold has zero volume in 
the original Hilbert space, the entropy, being the loga- 
rithm of this volume, diverges. 




FIG. 2: (ttab) as a function of the inverse temperature. No- 
tice the value {ttab) = 2/7V at j3 = (typical states). In the 
j3 — > oo limit we find the minimum (ttab) = 1/N. The phase 
transitions described in the text are at /3_ = —2/27, (tvab) = 
9/4N (left point) and f3+ = 2, (ty A b) = 5/4ZV (right point). 

With the same techniques, starting from (|3| we can 
find the cumulants of the purity for unbalanced biparti- 
tions. Leaving the details for a forthcoming publication 
we report the results for the first five cumulants only: 



1 2 



v^ab) 



/' 



8 2 + fx 



1 



iV 4 (i + M) 2 ' 
48 6 + 6^ 



384 22 + 33/z + 13/x 2 + ^ 3 



iV 13 



(1 + M) S 



(32) 



where /i = (M — N)/N. For /i = these reduce to the 
results of the previous section. 

Conclusions. We have calculated the generating 
function of a typical entanglement measure, averaged 
over the Hilbert space. We have shown that, when in- 
terpreted as a partition function, it possesses multiple 
phase transitions. In the different phases the distribution 
of Schmidt coefficients have different profiles. Sudden 
changes of these profiles occur at the phase transitions. 

We have studied these phase transition(s) as a function 
of a fictitious temperature (3, introduced to define the 
generating function of the purity. This fictitious temper- 
ature can also be thought of as localizing the measure on 
set of states with entanglement larger or smaller than the 
typical one |12] (in the same way temperature is used in 
classical statistical mechanics to fix the energy to a given 
value in the thermodynamic limit). 

Notice that the phase transitions investigated here, 
that appear in the study of the generating functions of 
any entanglement measure, are not quantum phase tran- 
sitions (QPT). Since entanglement is known to be a good 
indicator of QPTs [13] , it would be interesting to investi- 
gate the link, if any, between these different transitions. 

In conclusion, by using techniques borrowed from the 
study of random matrix theory, we gave a complete char- 
acterization of the statistics of one entanglement mea- 
sure. We also proposed one direction in which random 
matrix theory is likely to play a significant role in the 
study of entanglement, namely the role of the phase tran- 
sitions found in random matrix theory as describing the 
change in the profile of typical, less or more entangled 
states. 
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